Wave-number Selection by Target Patterns and Side Walls in Rayleigh-Benard 

Convection 



John R. Royer, Patrick O'Neill, Nathan Becker, and Guenter Ahlers 
Department of Physics and iQUEST, 
University of California, Santa Barbara, CA 93106 
(Dated: February 8, 2008) 

We present experimental results for patterns of Rayleigh-Benard convection in a cylindrical con- 
tainer with static side-wall forcing. The fluid used was methanol, with a Prandlt number <x = 7.17, 
and the aspect ratio was F = R/d ~ 19 (R is the radius and d the thickness of the fluid layer). 

In the presence of a small heat input along the side wall, a sudden jump of the temperature differ- 
ence AT from below to slightly above a critical value AT C produced a stable pattern of concentric 
rolls (a target pattern) with the central roll (the umbilicus) at the center of the cell. A quasi-static 
increase of e = AT/AT C — 1 beyond ei jC — 0.8 caused the umbilicus of the pattern to move off 
center. As observed by others, a further quasi-static increase of e up to e = 15.6 caused a sequence 
of transitions at Si t b, i = 1, 8, each associated with the loss of one convection roll at the umbilicus. 
Each loss of a roll was preceded by the displacement of the umbilicus away from the center of the 
cell. After each transition the umbilicus moved back toward but never quite reached the center. 
With decreasing e new rolls formed at the umbilicus when e was reduced below Ei t a < Ei,h- When 
decreasing e, large umbilicus displacements did not occur. 

In addition to quantitative measurements of the umbilicus displacement, we determined and 
analyzed the entire wave-director field of each image. The wave numbers varied in the axial direction, 
with minima at the umbilicus and at the cell wall and a maximum at a radial position close to 2F/3. 
The wave numbers at the maximum showed hysteretic jumps at e^j, and Ei, a , but on average agreed 
well with the theoretical predictions for the wave numbers selected in the far field of an infinitely 
extended target pattern. To our knowledge there is as yet no prediction for the wave number selected 
by the umbilicus itself, or by the cell wall of the finite experimental system. 

PACS numbers: 47.54.+r,47.27.Te,47.20.-k,47.20.Bp 



I. INTRODUCTION 

When a temperature difference AT exceeding a critical 
value AT C is applied across a thin, horizontal layer of 
fluid, convection occurs. This phenomenon is know as 
Rayleigh-Benard convection (RBC). The fluid flow then 
forms a pattern. The onset of convection is determined 
by the Rayleigh number 

R = agd^AT/nv . (1) 

Here a is the isobaric thermal expansion coefficient, g 
the acceptation of gravity, d the thickness of the sample, 
k the thermal diffusivity, and v the kinematic viscos- 
ity. Linear stability-analysis for a laterally infinite sys- 
tem shows that convection will occur with a critical wave 
number k c = 2tt / X c = 3.117 (A is measured in units of d) 
when R exceeds R c — 1707.3. For the laterally infi- 
nite system the patterns that evolve beyond the onset of 
convection and their stability are determined by R and 
by the Prandtl number 

o = v/n . (2) 

Weakly nonlinear theory predicts that immediately above 
onset, a laterally infinite uniform system should yield a 
pattern consisting of parallel, straight rolls. [2| For real 
experimental systems the patterns are influenced by the 
lateral boundaries which contain the fluid. Even in these 
finite systems excellent approximations of the predicted 



ideal straight rolls of the infinite system can be found 
under some conditions. 0,0,0 However, the boundaries 
can also lead to patterns of different symmetry. Thus, 
by heating a thin cylindrical convection cell gently from 
the side in addition to heating primarily from below, a 
pattern of concentric rolls can be stabilized. A concentric 
pattern can also be the result of horizontal temperature 
gradients near a cylindrical side wall which are intrinsic 
to the particular construction of the sample cell. The 
concentric pattern will consist of n rolls, where n is an 
integer. Of these, the one in the middle really would be 
more properly viewed as a convection cell, with either 
up-flow or down-flow in its center. It is often referred to 
as the umbilicus. This type of pattern, know also as a 
target pattern, was studied extensively by Koschmieder 
and Pallas [|| for relatively large a. The rolls have a 
mean dimensionless wave length A = 2T/n where T = 
R/d (R is the sample radius) is the aspect ratio of the 
sample and where n is the number of rolls along a radius. 
The corresponding wave number is k = 2ir/\ = irn/T. 
Koschmieder and Pallas found that the wave length of 
the rolls away from the center and the side wall increased 
with increasing AT. They also observed a sequence of 
transitions with increasing AT, each of which involved 
the loss of a roll at the sample center. 

Later work for smaller a showed that at a certain value 
£i jC of e = R/Rr. — 1 the target pattern becomes unsta- 
ble. Hi, H llfl Jnj The instability is known as the 
focus instability. 12] How this occurs in detail depends 
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on a, r, and on the nature of the confining side walls. 
Here we describe what happens for our particular sam- 
ple (u = 7.17, r = 19), where the side walls fixed the 
phase of the pattern at the radial position r = T (r is 
scaled by the cell thickness d). For the work of Ref. 8] 
((7 = 6.1, r = 7.5) the observed phenomena were simi- 
lar. As e is increased in small steps from small values, 
the umbilicus shrinks. The result is a gradual decrease 
of the local wave numbers of the rolls away from the 
center. When e is increased beyond ei !C , the umbilicus 
breaks the cylindrical symmetry by moving off center. 
This transition is continuous and the off-center patterns 
are stationary in time, but the radial umbilicus position 
is e dependent. As e is further increased, the umbilicus 
becomes even smaller, and then disappears at e = £\,b- 
At that point the pattern once more approaches a state of 
cylindrically symmetric rolls, although perfect rotational 
invariance is not fully recovered. With further increase 
of £, this process repeats itself, with another roll loss at 
£2,b > £i, 6 ; and so forth. The net result is a gradual re- 
duction with increasing e of the wave numbers of the rolls 
between transitions and a discontinuous change of k at 
each transition. The transitions are hysteretic, occurring 
at £i ia < £i,b when AT is decreased. 

We studied target patterns for < e < 15.6 in a sam- 
ple with r — 19.0 using methanol with a — 7.17. Wc 
developed high-resolution numerical umbilicus-detection 
algorithms and determined quantitatively the umbilicus 
displacement 6 with increasing as well as decreasing e. 
At small £ the results revealed the location of the focus 
instability, at e ~ 0.8. Previous measurements, [ll| for 
a = 0.93 and a larger T = 43, had found the instability at 
lower values, near e ~ 0.1. Over the whole e range of our 
experiment the data for 5, together with measurements 
of an average wave number (k) away from the center and 
the side wall, revealed eight hysteretic transitions. We 
observed large umbilicus displacement only for increasing 
e; for decreasing e the hysteresis yielded transition points 
£j iQ which were sufficiently low to avoid the £-range over 
which the umbilicus displacement was large. 

For £ > 5.6 cross rolls [HEi formed at the outermost 
roll when the system was close to a transition at e^. 
This was particularly pronounced on the side opposite to 
the umbilicus displacement where the local wave lengths 
were exceptionally large when 5 was large. Nonetheless, 
the general nature of the pattern was maintained up to 
our largest e values. For e > 14 periodic time depen- 
dence associated with the oscillatory instability [13, [l4j 
developed in the region near the cross rolls. 

A complete characterization of the pattern involves a 
knowledge of the entire wave-director field k(r, 9) where 
9 is the angular and r the radial position in the sample. 
Rather than simply measuring an average wave number 
along a cell diameter (as was done in previous work), we 
implemented a local wave-director analysis [l5l [To . Il7j 
and determined k(r,9). From k(r,9) we could then de- 
termine various averaged quantities, including (A;). The 
results for (k), determined over a radius range that ex- 



cluded the rolls near the center and the side wall, were 
discontinuous at the transitions. The hysteresis loop of 
each transition was traced out quantitatively. We found 
that the azimuthal average kg(r) = (k)g had interesting 
structure as a function of r, showing different selection 
at the umbilicus, in the bulk of the sample, and at the 
cell wall. The competition between the selection mecha- 
nisms yielded a broad maximum of kg(r) near r = 2T/3. 
For our experimental conditions it did not yield any ra- 
dially traveling waves, as predicted by Tuckerman and 
Barkley [^j for the case of conducting side walls. 
As was found in previous work, 0, H, EJ Ej the values 
of (k) selected in the bulk were in good agreement with 
the predicted wave numbers for the far field of infinitely 
extended concentric rolls when the discontinuities at the 
transitions were smoothed out. We also report results for 
the wave numbers selected at the umbilicus and at the 
side wall; but for these there seem to be no theoretical 
prediction. 

In the next section we shall review theoretical predic- 
tions for the selection by concentric rolls and relevant 
previous experiments. Then, in Sect. II I II we shall dis- 
cuss the experimental apparatus and procedures, as well 
as the image-analysis methods. Section 17x71 gives our re- 
sults. It consists of a discussion of the patterns observed 
in various £ ranges, of a presentation of the umbilicus- 
displacements results, and of a presentation of our wave- 
number results. A brief section summarizing this work 
ends the paper. 



II. THEORETICAL PREDICTIONS AND 
PREVIOUS EXPERIMENT 

Theoretical predictions of the wave numbers fcs(e) se- 
lected in the far field of infinitely extended target p at- 
terns have been made by several authors. [12 . 12C1 l2ll 122 . 

For the infinite system wave-number adjustment can 
take place by expansion of the pattern to large distances. 
Finite laboratory systems differ from this in an impor- 
tant way. The side wall, under typical experimental con- 
ditions, pins the phase of the pattern and prevents this 
expansion. Thus, a change of the average wave number 
with changing e is possible only at the umbilicus. How- 
ever, also at the umbilicus unhindered phase slip is not 
possible. Instead, a discontinuous and hysteretic process 
involving the destruction or creation of a convection cell 
occurs and leads to a discontinuous change of the wave- 
number field at the transition. The discontinuous effect 
on the average wave number decreases as the aspect ra- 
tio of the sample increases because the loss of a single 
cell in the center is a smaller perturbation for a larger 
sample. Thus, even for the sample with a boundary one 
expects, in the large-r limit, a continuous curve /cs(e) as 
a function of e. Theoretically the location of this curve 
is determined by the rotational symmetry of the target 
pattern which does not permit the mean flow that is in- 
duced by roll curvature under less symmetric conditions. 
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FIG. 1: The theoretical prediction for Sb = (fcs/fc c — l)/s 
in the limit of small e as a function of the Prandtl number 
a. The circle is the result of the resent work to be discussed 
below. 

Thus, the mean flow must be balanced precisely by a 
pressure gradient. [20l l2ll W\ This condition leads to a 
unique wave number ks(s). Initial predictions of the se- 
lected wave numbers were applicable only for small e and 
are given bv plE^ 

k B /k c -l = S B e + 0{e 2 ) (3) 

with S B = -N'/R 2 , N' = 0.1659 + 1.426/er - 1.220/cr 2 , 
and R 2 = 10.76 - 0.073/cr + 0.128/er 2 e. Here k c = 3.117 
is the critical wave number at the onset of convection. In 
Fig. \I\we show Sb as a function of a. For the present 
work we have a = 7.17 and Sb = —0.0317. The mea- 
surements to be described below yielded the experimental 
value —0.0285 ± 0.001. This result is shown as the circle 
in the figure. It falls slightly above, but is generally in 
good agreement with, the prediction. For the work of 
Ref. = 14) the theory yields S B = -0.0243, but no 
corresponding result has been extracted from the data. 
Measurements for a = 0.93 were reported in Ref. ^JJ an d 
yielded a value for Sb close to zero, but a precise number 
was not quoted. For that case one obtains the prediction 
Sb = —0.0266, but «Sb((7) changes rapidly with a and 
passes through zero at a — 0.79. 

Calculations valid at larger e were made by Buell and 
Catton [23| (BC) and by Newell et al. 0] using a com- 
bination of analytic and numerical methods. BC gave 
results for a few values of a and R, and we shall com- 
pare our measurements with interpolations between these 
predictions. We show their results for a Prandtl num- 
ber close to that of our experiment in Fig. [21 Also 
shown in that figure are the nearby stability bound- 
aries of infinitely extended uniform straight rolls. 0, 0] 
Although in the large-cr limit it is expected that the 
zig-zag instability will coincide with the selected wave 
number. [20I . l2ll |22| it is clear from the figure that a = 7 




FIG. 2: Stability boundaries of infinitely extended straight 
rolls for a Prandtl number a — 7.0. Dashed line: zigzag 
instability jl4^ . Solid line: cross-roll instability [Til ]. Also 
shown, as a dotted line, are the wave numbers predicted for 
the far field of infinitely extended concentric rolls [2!! . The 
two solid circles correspond to our experimental observation 
of the cross-roll instability near the side wall. 



is still far from that limit. 

To our knowledge there is as yet no prediction of the 
wave-number selection at the umbilicus itself. The results 
to be presented here show that the wave number selected 
there is smaller than the one in the far field. 

Aside from the phase pinning mentioned above, the 
side walls of experimental cells have an additional influ- 
ence on the pattern. The boundary conditions imposed 
by the walls can yield a separate and competing selec- 
tion mechanism. [pa, lifil |24| This mechanism is expected 
to depend on the conductivity of the side walls relative 
to that of the fluid, and it is not understood quantita- 
tively for typical experimental conditions. For perfectly 
conducting side walls, Tuckermann and Barkley [HI fl9| 
predicted a pattern of radially-traveling waves. The trav- 
eling nature of the pattern can be understood in terms of 
the competition between the selection by the curved rolls 
on the one hand and by the wall on the other which can 
lead to a wave-number gradient which in turn can lead 
to a non-zero time derivative of the phase of the pattern. 
[25L |26| However, to our knowledge such radially trav- 
elin g w aves have not yet been observed in experiments. 

mil 

For a values not too small the wave numbers predicted 
for curved rolls do not depend very strongly on a. In this 
a range the wave numbers of target patterns have been 
determined before, 0, 0> S @ but to our knowledge they 
were never compared quantitatively with the predictions 
of BC. Koschmieder and Pallas 6] determined a weighted 
average of the wave lengths of all rolls except for the 
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outermost one. In Fig. |2K we show their results for T = 
13.28, and for a — 511 and 916. The solid line is the 
prediction of Ref. |23j for a — oc. Aside from a lateral 
shift, the overall shape of a curve passing through the 
data is very similar to the theoretical curve. Since the 
data do not extrapolate to k c = 3.117 as e vanishes, the 
shift presumably is due to experimental uncertainties of 
the length scales involved in the determination of A. 

In Fig. |3) we show results from Ref. Q which were for 
a = 6.1 and T = 7.5. These data were obtained by mea- 
suring the average wave length of 2 rolls near r ~ 2r/3. 
They seem to show considerable scatter, but actually this 
is due to the discreteness of the values of k which is no- 
ticeable in the figure and which becomes more apparent 
for this relatively small T where the system contains only 
a small number of rolls. In this system hysteretic tran- 
sitions were clearly observed. The theoretical curve j2^| 
for this a value (solid line) is a good smoothed represen- 
tation of the data. 

Finally, in Fig. 02;, some of the results reported by Cro- 
quette and Pocheau [?J are shown. Those measurements 
were for a = 14 and T — 20, and were obtained by mea- 
suring the average roll width along a diameter in a region 
away from the center and sides where the roll width ap- 
peared constant. The solid line is the numerical result 
[23| for a = 14. The agreement clearly is very satisfy- 
ing. The data do not reveal much difference when e is 
increased (open symbols) or decreased (solid symbols). 
The authors report observing hysteretic transitions, and 
thus it is somewhat surprising that this hysteresis does 
not manifest itself in the wave-number measurements. 

All of the previous investigations revealed that the roll 
adjacent to the side wall was exceptionally wide, but none 
of the prior investigations made any attempt to deter- 
mine quantitatively the wave number selected at the wall 
or at the cell center. 
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III. EXPERIMENT 
A. Apparatus 

We used a Rayleigh-Benard convection apparatus de- 
scribed in detail elsewhere [2^, |3(|. The top plate 
was a single-crystal-sapphire disk, and the bottom plate 
was a polished aluminum disk with a metal-film heater 
mounted under it to provide uniform bottom-plate heat- 
ing. Temperature-controlled water was circulated over 
the top plate. The bottom-plate temperature was con- 
trolled to create the desired temperature difference. The 
mean temperature was maintained at 22.0°C. For both 
the top and bottom plate, the temperature varied about 
the set temperature by less then 10~ 3 °C. 

The cell wall was made of Lexan (thermal conductivity 
0.23 W/m K) with an outer diameter of 9.6 cm and an 
inner diameter of D = 8.89 cm. A 0.012 cm diameter 
manganin wire was embedded in the cell wall to provide 
the side- wall heating. This wire had a resistance of about 
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FIG. 3: Previous wave- number measurements for patterns of 
concentric rolls, (a): from Ref. for a — 511 (circles) and 
916 (squares), and for F = 13.28. For (b) and (c) open (solid) 
symbols were taken with increasing (decreasing) e. (b): from 
Ref. @ for a = 6.1 and T = 7.5 . (c): from Ref. for a = 14 
and r = 20. The solid lines in all figures are the predictions 
from Refs. 12311 and [12] for the relevant a values. 



13 f2, and typically dissipated 0.26 W. 

The height of the side wall was 0.229 cm, though the 
actual fluid height d depended on the compression of an 
O-ring which sealed the cell. The actual fluid height was 
determined by measuring AT C and then inferring d from 
the fluid properties and the measured AT C = 0.873 ± 
0.005°C. We found d = 0.234cm, giving an aspect ratio 
r = 19.0 . 

The fluid used was methanol which, at 22.0°C, had a 
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FIG. 4: The central 11. Ad x 11. Ad section of the pattern for 
e values near the transition at £4,5. They are from a run in 
which AT was increased with 0.05° C steps (Ae = 0.057). The 
system was equilibrated for one hour after each temperature 
step before taking a picture. The images were divided by a 
background image, and then processed numerically so as to 
greatly enhance their contrast. They are for (a) e = 3.139, 
(b) £ = 3.244, (c) £ = 3.301, and (d) e = 3.359. 



Prandlt number a = 7.17 and a thermal conductivity of 
0.20 W/m K. The vertical diffusion time was t v = 54s. 

The convection rolls were imaged using the shadow- 
graph method [^112 • At large e, strong spatial varia- 
tion of the refractive index caused nonlinear optical ef- 
fects in the images, [32j but the general features of the 
pattern could still be discerned. This effect limited our 
ability to do wave-number measurements for e > 4. 
Nonetheless, the location of the umbilicus could still be 
measured. 

In order to obtain concentric rolls, AT was set to zero 
and the side-wall heater-power was set to 0.26 W. After 
waiting two hours, we increased the temperature differ- 
ence to AT = 0.78°C (e ~ -0.1) and allowed the sys- 
tem to equilibrate for one hour. After this the temper- 
ature difference was abruptly increased to AT = 1.3° C 
(e ~ 0.5) and the system was allowed to equilibrate for 
two hours. The jump of AT was necessary because de- 
fects typically formed in the interior when the tempera- 
ture difference was raised gradually. Once concentric rolls 
were obtained, the temperature difference was adjusted 
in small steps (typically between 0.05°C and 0.01°C) and 
the system equilibrated for one hour before taking an im- 
age at each temperature step. The power supplied by the 
side-wall heater was kept constant throughout the run. 



B. Image Analysis 

Shadowgraph images of some of the patterns are shown 
in Figs.01and|3J All images were divided by a background 
image taken at AT = 0. They were filtered by setting 
to black (white) all pixels that fell below (above) a cer- 
tain threshold. A gaussian blur was then applied to the 
chopped image. 

We developed an algorithm to determine the local 
wave-director field k. This method will be described in 
detail elsewhere. First the orientation field 6{x, y) at each 
point of the image was calculated using a method similar 
to one introduced by Cross et. al., ^| but higher angular 
resolution was employed. The wave-number field k(x,y) 
was then calculated at each point by moving orthogonal 
to the the roll orientation at that point and determining 
the number of pixels which falls within one wave length, 
assuming a locally periodic structure. 

Determinations of the wave number 2ttcI/X from the 
images require a knowledge of d and of the horizontal 
distance Sx between adjacent pixels. The uncertainty of 
d yielded an uncertainty of only a small fraction of a 
percent. We determined Sx by counting the number of 
pixels which spanned the sample diameter. This yielded 
Sx = 0.095 ± .002 and led to a combined uncertainty of 
about 2 % for k. An extrapolation of the measured k- 
values with decreasing e to e — gave k c = 3.166, which 
is within the experimental uncertainty of the value 3.117 
for the infinitely extended set of straight rolls. 

We measured the location of the umbilicus by using a 
24d x 2Ad square section surrounding the center. The lo- 
cation of the umbilicus was estimated by computing the 
maximum of the cross-correlation of the measured orien- 
tation field 0(x,y) and the orientation field for perfectly 
centered rolls 9 c (x,y). This estimate was used only to 
identify the rough location of the umbilicus. Once the 
umbilicus was identified, we used an edge-detection algo- 
rithm on the umbilicus to determine its position with a 
resolution of about 70^m or 0.03d. 

IV. RESULTS 

A. Patterns: Increasing e 

We covered the range < e < 15.6 and with increas- 
ing e observed eight transitions at e^i = 1, ...,8 under 
quasi-static conditions. A closeup view of images near 
the third transition is shown in Fig. At each transi- 
tion a roll was lost at the umbilicus. For all but the first 
transition the umbilicus first moved away from the cen- 
ter as e grew. The angle of the umbilicus displacement 
did not have a preferred value, as can be seen from the 
examples in Figs. 0] and |SJ This indicates that the cell 
was sufficiently uniform. At the transition the umbilicus 
collapsed, reducing the number of rolls in the system by 
one. After a roll was lost, the new umbilicus immediately 
moved back toward the center of the sample but never 
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FIG. 5: Images near the transitions at £6,6 and £7^. The 
optical effects that distort the roll patterns can be seen in 
images (b) and (e) in the rolls that are being compressed by 
the off-center umbilicus. These images are from the same run 
as those in Fig. El The values of e are (a) 5.649, (b) 5.707, 
(c) 5.821, (d) 6.910, (e) 7.368, and (f) 7.482. 



quite reached it. This process was analogous to that ob- 
served in Refs. [7j and |8|]. The umbilicus displacement 
before each transition produced a noticeable azimuthal 
variation of the wave number, with higher wave num- 
bers in the direction of the displacement, and lower wave 
numbers opposite the displacement direction. 

Close to but before the sixth transition, for e ~ 5.4, a 
patch of cross rolls [lj| formed along the outermost 
roll on the side opposite the direction of umbilicus dis- 
placement. At that point the local wave number k of 
the concentric rolls was the smallest of the entire wave- 
number field, and was about 1.95 (k/k c — 1 = —0.39). 
The corresponding point in the e — k plane is shown 
as a solid circle in Fig. [21 The wave number b of the 
cross rolls that formed was about 4.3 . The value of k 
is somewhat larger than the prediction ^4| kca = 1-67 
(kcnjkc — 1 = —0.46) for laterally infinite straight rolls 
at this value of e. This suggests that the roll curvature 
or the side wall reduces the stablity of the rolls against 
the cross-roll perturbation. The wave number b of the 
cross rolls that formed is larger than the predicted value 
[pjj ben — 3.5. The cross rolls disappeared when e was 
increased further, the pattern lost the middle roll, and 
the umbilicus returned toward the center. At that point 




FIG. 6: Images for (a) e = 12.75 and (b) e = 15.56 from the 
same run as the one used for Figs. 0]and|3. The pattern in 
image (a) was stationary, while in image (b) the cross rolls 
oscillated along the axis of the main rolls. At these large e 
values the images are strongly influenced by nonlinear effects 
in the shadowgraph method. 



all parts of the pattern had returned to wave numbers 
safely in the stable Busse balloon. Before the next (sev- 
enth) transition, near e ~ 6.9, cross rolls appeared again. 
In this case we found k =1.86 (k/k c — 1 = —0.41) and 
b = 4.3. This point is shown as well in Fig. [21 Again the 
value of b is somewhat larger than the prediction b ~ 3.9 
for infinitely extended uniform straight rolls. After this 
transition the cross rolls fanned out to cover the entire 
outer roll. Images near the two transitions are shown in 
Fig. 

The cross rolls remained along the outer roll, and after 
the eighth transition, spread to the inner rolls as e was 
further increased. Despite the cross rolls, there was still 
a discernible pattern of nearly concentric rolls, with the 
cross rolls superimposed on these. Around s ~ 11.5 the 
focus began to move off center but did not complete a 
ninth transition. This is illustrated in Fig. [SJa for e = 
12.75. Cross rolls spread to the inner rolls, and as e was 
further increased, the umbilicus returned to the center 
without losing a roll, as seen in Fig. 03 for e = 15.56. 

The pattern became time dependent near e = 14.2, 
where small patches of traveling waves developed. These 
were superimposed upon and orthogonal to the concen- 
tric roll, and traveled along the axes of the concentric 
rolls. The travelling waves occurred first on the third or 
fourth roll from the wall, and with increasing e spread 
throughout the cell. We asociate this phenomenon with 
the oscillatory instability predicted by Clever and Busse, 
[l3L Il4| although for our Prandtl number the cross-roll 
instability precedes the oscillatory instability. 

The cross rolls observed here had been previously seen 
by Croquette et. al. Q, though there are differences in 
the way they appeared. Croquette et. al. observed that 
cross rolls moved in from the side wall only after e became 
larger than about 10. They did not observe a relation- 
ship between the umbilicus transitions and the cross-roll 
formation. Despite nearly identical Prandlt numbers and 
aspect ratios {a = 7, T = 20), Croquette et. al. could 
produce off-center patterns near transitions comparable 
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FIG. 7: The central 11. 4d x 11. 4d sections of images obtained 
with decreasing e near the transition at E4. a - This is from a 
run where e was first increased to 3.49 and then decreased 
in steps Ae = -0.01. The values of e are (a) 3.49, (b) 2.94, 
(c) 2.67, and (d) 2.65. Note how the size of the middle roll 
changes from image (a) to (c), and how in (c) the umbilicus 
remains close to the center even for an e value only just above 

£4,a- 



to those shown in Fig. [S] with no cross rolls In our 
system cross rolls permanently covered the outer roll for 
e = 7.48, well before Croquette et. al. first observed any. 
This difference can perhaps be explained by a difference 
in the strength of side-wall forcing. A future study of 
the influence of various levels of side-wall heating should 
shed some light on this issue. 



B. Patterns: Decreasing e 

The nature of the patterns for decreasing e depended 
on the initial state. When starting above the oscillatory 
instability (say e > 14), decreasing e broke up the concen- 
tric rolls. Wall foci evolved, displaced the umbilicus, and 
emitted new rolls as e was further decreased. However, 
the side-wall forcing was strong enough to re-orient the 
outermost rolls parallel to the sidewall for e <■ 0.21. By 
e ~ 0.09 the outer roll had returned to being completely 
parallel to the sidewall. However, defects remained in the 
interior of the pattern, and the system did not return to 
concentric rolls as e was decreased below zero. 

Starting from e < 6 with concentric time-independent 
rolls, new rolls were generated at the umbilicus for dis- 
tinct values £i^ a < Ei.b as e was decreased quasi-statically. 
During these transitions the umbilicus stayed close to the 
center, as illustrated in Fig. [7] We did not determine the 
maximum e at which we could start decreasing e without 
loosing the concentric pattern. 



FIG. 8: The umbilicus displacement S away from the sample 
center as a function of e for increasing e. This data is from 
the same run as the images in Figs. 0] and |S] where we in- 
creased e in steps of Ae = 0.057 and waited 1 hour after each 
temperature step for the system to equilibrate before taking 
a picture. 

C. Umbilicus Displacement 

With increasing e we measured the displacement S of 
the umbilicus from the center of the cell as a function 
of e for 0.5 < e < 12. Results over the entire e-range 
are shown in Fig. |SJ A more detailed plot of the re- 
sults for e < 4 is shown in Fig. |5| Seven transitions, at 
Ei b,i = 2,..., 8, are apparent from these data. At the 
first transition, at E\ t>i we did not resolve any discon- 
tinuity in 5. Prior to each observable transition there 
was a large displacement, followed by a relaxation back 
toward a centro-symmetric pattern when the umbilicus 
collapsed. 

Figure ITUI shows the data for 5 over the range e < 2. 
Here we see that the pattern is, within possible system- 
atic errors of 8, accurately centered for e < £x,c ~ 0.8. 
Above £i jC the increase of 8 is consistent with an initially 
linear dependence, and thus we fitted the results to 

(5(e) = S +di(e - £i, c ) + S 2 (e - e hc ) 2 (4) 

over the range £i jC < e < 1.7 and to 5(e) — So for e < 
£i iC . Here So corresponds to the small offset, well within 
our possible systematic errors, which is found even at 
small epsilon. The parameters 5o,5\,S2, and £i jC were 
least-squares adjusted. The fit gave £i. c = 0.86 ± 0.13, 
5i = 0.077 ± 0.075, and 5 2 = 0.12 ± 0.05. It is shown 
as the dashed curve in Fig. ^| The statistical error of 
S\ indicates that a fit to a quadratic equation (i.e. Eq.Q] 
with Si = 0) should be equally good. It yielded £i jC = 
0.66 ± 0.08 and S 2 = 0.14 ± 0.02. We identify e hc as the 
focus instability, i.e. as the first instability of the centro- 
symmetric pattern with our prevailing wave number and 
our Prandtl number a = 7.17. 

When £ was decreased, transitions involving the addi- 
tion of a roll at the umbilicus occurred at £i jQ , but were 
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FIG. 9: An expanded view at relatively small e of the results 
shown in Fig. |H] The thin solid line shows the umbilicus 
displacement for decreasing e. 




FIG. 10: A very expanded view at small e of the results shown 
in Figs. |H| and [5] This graph reveals the initial focus instabil- 
ity of the centro-symmetric pattern at ei jC . The points ei, a 
and £1,6 are the limits of the hysteresis loop associated with 
the transition which are revealed by the wave-number mea- 
surements to be discussed below. The dashed curve is a fit of 
Eq.0]to the data. 



not associated with large displacements of the umbilicus. 
Experimental results for S with decreasing e are shown by 
the thin solid line in Fig. This difference between in- 
creasingand decreasing e was observed also by Steinberg 
et. al. [|. 



D. Wave-Number Measurements 

We measured the averaged wave numbers (k) for < 
e < 4 for both increasing and decreasing s. For compar- 
ison with theoretical predictions one would like to make 
this measurement well away from the umbilicus; in prac- 
tice there is a limit set by the aspect ratio of the sample 
and by a competing selection mechanism associated with 
the side wall. As we shall discuss below, the local wave 




FIG. 11: A chopped and blurred image obtained with increas- 
ing e at e — 2.61. The highlighted section covers the range 
10 < r < 14 and shows the annulus over which we averaged 
the wave numbers to obtain the values (k) shown in Figs. 1121 
[T31 and ITU below. 




FIG. 12: Results for (k) obtained with decreasing e. The solid 
line is a fit of Eq.|K|to the data. 

numbers have a maximum at a radial position in the 
range 10 < r < 14 (0.53r < r < 0.74T). Thus we se- 
lected an annular region extending over this radial range 
as illustrated by the highlighted region in Fig.^Jto com- 
pute (k). In Fig. E| the results obtained with decreasing 
and relatively small e are shown. The solid line is a fit of 
the equation 

(k)/k c -l = S B e + S 2 e 2 (5) 

to the data over the range e < 0.6. The fit yielded 
k c = 3.167, Sb = -0.0285 ± 0.0010, and S 2 = -0.027. 
The result for k c , within the a priori estimate of the 
experimental uncertainties, agrees with the theoretical 
value k c = 3.117 for infinitely extended straight rolls. 
The result for Sb is shown in Fig. ^ as an open circle. It 
falls slightly above the predicted value —0.0317. 




FIG. 13: Circles: the average selected reduced wave num- 
bers (k)/k c — 1 vs. e. Squares: the reduced wave numbers 
ke(T)/k c — 1 selected at the side wall. Open symbols: in- 
creasing e. Closed symbols: decreasing e. Solid line thr ough 
the data for (k)/k c — 1: the prediction of Refs. |23l| and (l2|. 
Lower solid line: the zigzag instability for a = 7. 

Results for (k) /k c — 1 over the entire e range are shown 
in Fig. E| for both increasing (open circles) and decreas- 
ing (solid circles) s. Within experimental uncertainties 
they are consistent with previous measurements for sim- 
ilar er. The jumps in the wave number correspond 
to the loss or creation of a roll at a transition. The hys- 
teresis in the wave-number selection, previously observed 
by others, 0,13 is clearly visible in our results. R indi- 
cates partial pinning of the phase of the pattern at the 
umbilicus. We believe that this phase pinning is responsi- 
ble also for the difference between experiment and theory 
for the value of Sb ■ On average the data are in excellent 
agreement with the calculations by Buell and Catton j^l 
which are given by the solid line through the data. 

In order to illustrate the relationship between the first 
hysteretic transition at S\, a and £i,t> and the focus insta- 
bility at £i, c , we give an expanded view of the selected 
wave numbers at small e in Fig. One sees that the 
discontinuities at S\^ a and e\,b are quite sharp, and that 
£i jC is located near the middle of this first hysteresis loop. 
Values of Ei tll and Si^ for all eight transitions are given 
in Table □ 

The top part of Fig. 1151 gives an example for the radial 
variation of the shadowgraph intensity. In the bottom 
part we display the radial variation of the azimuthal aver- 
age kg of k for three e values. Also shown for comparison 
(bottom curve, down-shifted by 0.1) is the result of our 
analysis for a synthetic pattern of concentric rolls with a 
radial wave-number distribution shown by the solid line 
which represents a smooth curve through the results for 



FIG. 14: A detailed view of selected wave numbers (k) near 
the focus instability at E\, c and the first hysteretic transition 
at £i :( , and ei, - The symbols are as in Fig. 1131 



TABLE I: Location of the hysteretic transitions with decreas- 
ing (£i, a ) and increasing (e^) e. For Ei, a , i = 1, .., 4 the 
uncertainty Se = 0.006. For £5, a and ee.a, Se = 0.11. For Ei t b, 
Se = 0.029. 

1 0.598 1.154 

2 1.217 1.841 

3 1.858 2.528 

4 2.660 3.330 

5 3.6 4.437 

6 4.8 5.759 

7 - 7.368 

8 - 9.552 



£ = 2.614. 

Our analysis yields results for kg(r) which oscillate as a 
function of r. This is particularly noticeable for r below 
the maximum of kg(r). These oscillations are not present 
in the analysis of a synthetic concentric pattern with a 
uniform wave number, but as shown by the lowest set in 
Fig. UU they do appear in the analysis of a concentric 
pattern with a radial wave-number distribution equal to 
the smooth curve through the data for e — 2.614. Thus 
we regard them to be an artifact of the numerical pro- 
cedures. We note, however, that the amplitude of the 
oscillations is only about 0.5% of fc, and that a smooth 
curve through them is a good representation of the actual 
wave numbers. 

From Fig. ^| one sees that the experimental data for 
kg, as mentioned above, have a maximum in the range 
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radius (d) 

FIG. 15: Top: an azimuthal average of the shadowgraph in- 
tensity (arbitrary scale) as a function of radial position for 
e = 2.61 . Bottom: azimuthal averages of the local wave 
numbers. The vertical dashed line shows the location of the 
side wall, and the two vertical dotted lines show the inner and 
outer edge of the annulus highlighted in Fig llll and used to 
compute the average wave numbers (k) shown in Figs. 1121 1131 
and 1141 From top to bottom, the first three data sets are for 
e = 0.077, 1.297, and 2.614. The lowest set is for a synthetic 
concentric pattern (down-shifted in the figure by 0.1) with a 
wave-number distribution (shown by the solid line) given by 
straight-line representations of the experiment for e = 2.614. 
The solid straight lines through the experimental data at large 
r are fits to the data with r > 14 and were used to find kg (F) . 
The dashed straight lines through the experimental data at 
small r are fits for 5 < r < 10 and were used to determine 
fc(0). 



10 < r < 14, and this range was used to compute (k). 
For r > 14 the data were fitted by straight lines as shown 
in the figure, and these fits were used to extrapolate the 
selected kg(r) to kg(r = T). The results for kg(T) are 
shown in Fig. ^5] above. They show, as noted by others, 
[HOHH that the outermost roll has an anomalously large 
wave length (small wave number) . We note that kg (T) is 
not influenced significantly by the transitions at the um- 
bilicus. Interestingly, the wall-selected wave number does 
not extrapolate to k c as e vanishes. Instead, it crosses 
the zigzag instability-line of the laterally infinite system; 
but in the experiment for this finite system no instabil- 
ity was observed. We are not aware of any quantitative 
previous measurements or predictions for kg(T). 

The experimental results at small r suggest that yet 
another unique wave number is selected by the umbilicus 
itself. Fits of straight lines to the data, this time with 
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k/k c - 1 

FIG. 16: Triangles: The reduced wave numbers k(0)/k c — 1 
selected at the origin. Open symbols: increasing e. Closed 
symbols: decreasing e. Rough solid line: kg(T) for decreasing 
e. Smooth solid lines: predictions for (k) |23| and for the 
zigzag and cross- roll instabilities [Til ]. 



5 < r < 10, could be used to estimate the selected values 
fc(0) at r = 0. The results are given in Fig. For 
comparison we show in that figure also the prediction 
for the concentric-roll selection [23 (smooth solid line) 
and the experimental results for kg(T) (jagged line). One 
sees that fc(0) is influenced strongly by the hysteretic 
transitions at the umbilicus. At large e fc(0) is on average 
close to the result for kg(T), but as e decreases, k(0) is 
somewhat larger and tends toward an intermediate value 
somewhere between kg(T) and (k) as e — » 0. Also for 
k(0) there appear to be no previous measurements or 
predictions. 



V. CONCLUSION 

We studied Rayleigh-Benard convection of a fluid with 
a Prandtl number a — 7.17 in a cylindrical container of 
aspect ratio T = 19 in the presence of static side- wall 
forcing. Patterns of concentric rolls were obtained, and 
we studied the quasi-static evolution of these patterns for 
both increasing and decreasing e. 

For increasing e, over a range 0.4 < e < 15.6, the 
pattern underwent eight transitions where the middle roll 
moved off center and then disappeared as e was increased 
quasi-statically. We measured the displacement of the 
umbilicus for 0.5 < e < 12. Above e ~ 14 the pattern 
became time dependent with oscillating rolls traveling 
along the axis of the concentric rolls. 

For decreasing e, the concentric-roll pattern was lost 
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when the initial e was too large. When e was decreased 
from a moderate initial value, the concentric-roll pattern 
remained and new rolls nucleated at the umbilicus. Any 
umbilicus displacement preceding the nucleation of new 
rolls was modest, in sharp contrast with the behavior of 
the pattern for increasing e. 

Using new image- analysis techniques, we determined 
the wave-number field for < e < 4. From this we 
computed the azimuthal average kg. There was a very 
noticeable radial gradient of kg both near the side wall at 
r = r and near the center at r — 0. Averaging the wave- 
number field over the annulus 10 < r < 14 where kg was 
relatively constant, we obtained an average wave number 
(k). We found that (k) was discontinuous and hysteretic 
at the transitions where new rolls were formed or disap- 
peared at the umbilicus. Nonetheless, a smooth curve 
through the data agreed well with predictions [23] for 
the wave-number selection in the far field of concentric 



rolls. These average wave numbers were also consistent 
with previous work by others. 0, @ Near the wall and 
near the center the measured kg were extrapolated to 
r = r and to r = to obtain the wave numbers kg(T) 
and fc(0) selected by the wall and the umbilicus. We are 
not aware of previous measurements or predictions for 
these quantities. 
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